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PACS 87. 16.Dg - Membranes, bilayers, and vesicles 
PACS 87 . 15 . Vv - Diffusion 
PACS 68 . 35 . Np - Adhesion 

Abstract. - The adhesion of cells is mediated by membrane receptors that bind to complemen- 
tary ligands in apposing cell membranes. It is generally assumed that the lateral diffusion of 
mobile receptor-ligand bonds in membrane-membrane adhesion zones is slower than the diffusion 
of unbound receptors and ligands. We find that this slowing down is not only caused by the larger 
size of the bound receptor-ligand complexes, but also by thermal fluctuations of the membrane 
shape. We model two adhering membranes as elastic sheets pinned together by receptor-ligand 
bonds and study the diffusion of the bonds using Monte Carlo simulations. In our model, the 
fluctuations reduce the bond diffusion constant in planar membranes by a factor close to 2 in the 
biologically relevant regime of small bond concentrations. 
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Introduction. — Advances in optical microscopy have 
made it possible to observe the diffusion of single lipid 
and protein molecules in biological membranes with a po- 
sitional accuracy down to about 40 nm at ms time resolu- 
tion [1,2]. Deviations from free Brownian motion were fre- 
quently found [2,3] and interpreted as the consequence of 
membrane microheterogeneity. Moreover, cis-interactions 
of mobile membrane proteins with cytoskeletal-anchored 
proteins or aggregates were identified via transient im- 
mobilization [4]. Trans-interactions between single cellu- 
lar cadherin molecules and their soluble counterpart have 
been studied [5] , paving the way for a characterization of 
membrane-membrane adhesion zones at the level of indi- 
vidual adhesion molecules. 

Current single-molecule experiments aim at detecting 
the reduced diffusion constants of bound receptors and 
ligands as important parameters that report on the time 
course of interaction in membrane adhesion zones. It is 
generally expected that the diffusion of bound receptor- 
ligand complexes is slower than the diffusion of the con- 
stituent unbound receptor and ligand molecules, simply 
because the complexes are larger and therefore also ex- 
posed to a larger viscous drag. We focus here on an- 
other and, to our knowledge, so far neglected aspect of 
the lateral diffusion of receptor-ligand bonds: The cou- 
pling of the diffusion to thermal membrane shape fluctu- 



ations. We show that shape fluctuations on a nanometer 
scale significantly decrease the diffusion constant of the 
bonds. On this length scale, the fluctuations are governed 
by the bending energy of the membranes. 

We consider a statistical-mechanical model in which 
two apposing membranes are discretized into quadratic 
patches. Each of the patches in the first membrane can be 
occupied by a single receptor molecule, which can bind 
a ligand molecule in the apposing patch of the second 
membrane. We consider here irreversibly bound receptor- 
ligand bonds that impose a local membrane separation 
close to the bond length l . The lateral diffusion of the 
bonds is modeled as a hopping process to neighboring 
patches and implies that the local separation at these 
patches is changed to l a . Hopping steps therefore can lead 
to an increase of the membrane bending energy, which 
slows down diffusion. 

We show that the decrease of the diffusion constant de- 
pends on the thermal membrane roughness, which is a 
measure for the strength of fluctuations. The roughness, in 
turn, depends on the concentration of the receptor-ligand 
bonds: it decreases with increasing bond concentration. 
For small bond concentrations, the membrane roughness 
at neighboring membrane patches of receptor-ligand bonds 
tends towards a finite value. In our statistical-mechanical 
model, the fluctuations reduce the diffusion constant by a 
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factor of about 1.8 in the regime of small bond concen- 
trations, which is the biologically relevant regime for typi- 
cal concentrations of receptor-ligand bonds in membrane- 
membrane adhesion zones. 

Model. — We consider two on average parallel 
membranes that are interconnected by receptor-ligand 
bonds. In our model, the membranes are discretized into 
quadratic patches of size a x a. The linear patch size a 
corresponds to the smallest bending deformations, which 
are around 5 nm, according to simulations with molecular 
membrane models [7] . Each patch i in the first membrane 
directly apposes a patch in the second membrane. The 
membrane conformations then can be described by the lo- 
cal separation k of all pairs of apposing patches. Two 
apposing membrane patches can be pinned together by 
a single receptor-ligand bond. We assume here that the 
bonds are irreversible and rather rigid. At patch sites i 
with receptor-ligand bonds, the local separation is con- 
strained then to the length k = l a of the receptor-ligand 
bonds. 

The conformations of the two membranes are governed 
by the bending energy of the membranes, under the con- 
straints that the local separation is equal to l a at all 
patch sites i where receptor-ligand bonds are present. To 
simplify the notation, we use the rcscalcd separation field 
Zi = (lij 'a)y 'n/ksT , where ks is Boltzmann's constant, 
and T is temperature. Here, n — kik 2 /{ki + k 2 ) is the ef- 
fective bending rigidity of the two membranes with rigidi- 
ties K\ and K2- The conformations of the rescaled sepa- 
ration field Zi are governed by the effective bending en- 
ergy [8] 
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Fig. 1: (Left) A receptor-ligand complex that interconnects 
two planar and parallel membranes. The lateral diffusion of 
the receptor-ligand complex, indicated by the arrows, does not 
involve any membrane bending. - (Right) A receptor-ligand 
complex interconnecting two thermally fluctuating or 'rough' 
membranes. The lateral diffusion is coupled to the membrane 
fluctuations. 
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with A d Zi = za + z l2 + z i3 + z j4 - Azi where z lY to z iA 
are the membrane separations at the four nearest-neighbor 
patches of membrane patch i. The discretized Laplacian 
AdZi is proportional to the local mean curvature of the 
separation field. We focus here on the bending energy of 
the membranes since this energy is the dominant contri- 
bution to the elastic energy on the small length scales up 
to about 100 nm that are relevant here. Elastic contribu- 
tions from a membrane tension a or from the cytoskeleton 
anchored to cell membranes are only important on larger 
length scales, i.e. on length scales larger than the crossover 
length y/n/a, which has a typical size of several hundred 
nanometers for rigidities k between 10 and 20 fc^T [9] and 
tensions a of a few /j,3/m [10], and larger than the aver- 
age separation of neighboring cytoskeletal anchor points, 
which is about 100 nm [6]. 

The lateral diffusion of the receptor-ligand bonds is de- 
scribed as a hopping process on the square lattice of mem- 
brane patches. Let us first consider the hypothetical situ- 
ation in which the membranes are planar and have a con- 
stant local separation z; t = z Q where z Q = {l Q / a)^J n/hsT 
is the rescaled length of the receptor-ligand bonds. This 
situation is depicted in the left cartoon of Fig. 1. The lat- 



eral diffusion of the receptor-ligand bonds is then indepen- 
dent of the membrane conformations. During each time 
step td of the diffusive hopping process, each receptor- 
ligand bond can hop to one of the four nearest-neighbor 
patches, provided the patch is not occupied already by 
another receptor-ligand bond. At low concentrations of 
receptor-ligand bonds, diffusive encounters of two bonds 
can be neglected and the diffusion of a single receptor- 
ligand bond is characterized by a sequence of 'jump vec- 
tors' Gtj with length |eij| = a. The orientation of each 
jump vector along one of the four 'hopping directions' 
on the square lattice is independent of the orientations 
of the other jump vectors. After N jumps, the lateral 
mean-square displacement r 2 of r = a i then sim- 

ply is r 2 = Na 2 , and the diffusion constant D, defined by 
r 2 = 4Dt with t = Nt d , is D = a 2 /4t d . 

Lipid membranes, however, are not planar. They ex- 
hibit significant thermal shape fluctuations since their typ- 
ical bending rigidities are between 10 and 20 fc^T [9], 
i.e. only one order of magnitude larger than the thermal 
energy ksT. The fluctuations affect the lateral diffusion 
of the bound receptor-ligand complexes. The jumps of 
receptor-ligand bonds to neighboring membrane patches j 
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then depend on the local separation Zj of these patches. 
We assume here that a jump, or hopping step, to a neigh- 
boring patch j implies that the local separation Zj of this 
patch is shifted to z a , the rescaled length of the receptor- 
ligand bond. The hopping step then changes the bending 
energy (1) of the membranes. We study the diffusion in 
Monte Carlo simulations with the standard Metropolis cri- 
terion [11]. According to this criterion, a hopping step is 
always accepted if the bending energy is reduced by the 
step. If the bending energy is increased by AE, the step is 
accepted with probability exp[— AE/UbT]. Hopping steps 
are always rejected if the neighboring membrane patch is 
already occupied by another receptor-ligand bond. The 
thermal fluctuations of the membranes are simulated via 
Monte Carlo steps that shift the local separation Zi of 
patch i to a 'new' value Zi + p where p is a random number 
between -1 and 1. These local moves of the separation field 
are again accepted or rejected according to the Metropolis 
criterion, except for moves with z% + p < 0, which are al- 
ways rejected since the membranes cannot penetrate each 
other. 

An important aspect are the characteristic time scales 
for the lateral diffusion and for the membrane fluctua- 
tions on the smallest length scale a ~ 5 nm of the model. 
Typical diffusion constants D for membrane proteins are 
between 0.1 and 1 pm 2 /s [12, 13], which correspond to 
characteristic jump times td — a 2 /AD between 1 and a 
few microseconds. The characteristic relaxation time for 
membrane shape fluctuations on the length scale a ~ 5 nm 
is around 0.1 microseconds [14], i.e. one order of magni- 
tude smaller than the characteristic diffusion time on this 
length scale. To capture this difference in the two time 
scales, we perform on average 10 Monte Carlo steps for 
the local separation Zi on each lattice site between con- 
secutive hopping steps of the receptor ligand-bonds. The 
local separation at the neighboring patches of a receptor- 
ligand bond 'anneals' during these 10 Monte Carlo steps, 
i.e. the separations arc not, or only very weakly, time- 
correlated with the separations of these patches at the 
last hopping step of the receptor-ligand bond. 

The data presented here are from simulations with up to 
10 7 Monte Carlo steps per membrane patch for the sep- 
arations zi. We have simulated membranes with a size 
up to 120 x 120 patches and periodic boundary condi- 
tions. In our simulations, the lateral correlation length of 
the separation field is significantly smaller than the linear 
membrane size. Our results are therefore not affected by 
the finite membrane size. The initial conformation of the 
Monte Carlo simulations is the planar membrane with sep- 
aration Zi = z a and random distribution of receptor-ligand 
bonds. From this initial conformation, the membranes re- 
lax relatively fast into thermal equilibrium, for which the 
thermal averages reported below have been measured. 

Results. — Our model has only two parameters: (i) 
The area fraction X of receptor-ligand bonds, i.e. the frac- 
tion of membrane patches that contain bonds, and (ii) the 




Fig. 2: (a) Rejection probabilities for lateral diffusion steps, or 
'hopping steps', as a function of the area fraction X of receptor- 
ligand bonds in the limit of large rescaled bond length z . The 
dotted line represents the probability P ex that a step is re- 
jected because the membrane patch is already occupied by a 
receptor-ligand bond ('exclusion rejection'), the full line repre- 
sents the probability P/j that a step is rejected because of an 
increase in the membranes' elastic energy ('fluctuation rejec- 
tion'). The dashed line is the total rejection probability Ptot- 
The lines result from interpolation of the data points obtained 
in Monte Carlo simulations. The statistical errors for the the 
data points are smaller than the symbol sizes. - (b) Overall 
membrane roughness £± (dashed line) and roughness £x,n at 
unoccupied neighbor sites of receptor-ligand bonds (full line) as 
a function of the bond area fraction X. The lines are obtained 
from interpolation of the Monte Carlo data points. 



rescaled length z a = (l /a) ^n/ksT of the bonds, which 
depends on the actual bond length l D , the effective bend- 
ing rigidity n and the temperature T. The characteris- 
tic length scale related to the area fraction is the mean 
bond separation = l/VX- Other length scales, such 
as the average membrane separation (z^ and the mem- 
brane roughness £j_ = sj {zf} — (zi) 2 , depend on these two 
parameters (see below). Here, (...) denotes a thermal 
average over membrane conformations. 

We first consider the case where z D is significantly larger 
than the membrane roughness In the limiting case 
of large z a , the fluctuating membranes don't 'touch' each 
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Fig. 3: (a) The overall membrane roughness £x for large 
rescaled bond lengths z a versus the mean bond separation 
(,bb = 1/ \fX in units of the linear size a of membrane patches. 
Here, X is the area fraction of receptor-ligand bonds. The 
line, given by £x = 0.02 + 0.14^66, is obtained by a linear fit 
to the last four data points for large £w>. - (b) Fluctuation 
rejection probability Pfi for lateral hopping steps versus the 
roughness £x,n at unoccupied nearest neighbor sites of bonds 
on the quadratic array of membrane patches. At large values 
of £x,n, which correspond to small bond area fractions X, the 
rejection probability Pfi is proportional to £±,„. The dashed 
line Pfi — —0.24 + 1.35 £x, n is obtained by a linear fit to the 
four Monte Carlo data points for £x jn > 0.41. 



other, and the hard-wall repulsion at local separations 
Zi = docs not play any role. Therefore, in this limit 
the results do not depend on z a , i.e. the area fraction X 
of receptor-ligand bonds is the only parameter that affects 
the lateral diffusion of the bonds and the membrane rough- 
ness. In addition to the overall membrane roughness £± 
defined above, we will also consider the membrane rough- 
ness £x,n at unoccupied nearest-neighbor sites of bonds 
on the quadratic array of membrane patches. These un- 
occupied nearest-neighbor sites are the sites to which the 
bonds can jump during diffusion. 

The effect of the membrane fluctuations on the lateral 
diffusion of the receptor-ligand bonds can be character- 
ized by the rejection probabilities of hopping steps. There 
are two reasons why hopping steps can be rejected. First, 



a hopping step is rejected if the neighboring membrane 
patch is already occupied by a receptor-ligand bond. We 
call this 'exclusion rejection', since the bonds exclude each 
other from the patches. The second reason is the 'fluctu- 
ation rejection', i.e. hopping steps are rejected with prob- 
ability cxp[— AE/ksT] if the step increases the bending 
energy of the fluctuating membranes by AE. These two 
rejection probabilities and their sum, the total rejection 
probability, are shown in Fig. 2(a) as a function of the 
bond area fraction X in the limit of large z . The exclu- 
sion rejection (dotted line) increases with the bond area 
fraction. The fluctuation rejection (full line), in contrast, 
decreases with the area fraction. The total rejection prob- 
ability has a minimum at intermediate area fractions. 

Concentrations of receptor-ligand bonds in cell-cell ad- 
hesion zones can be as large as several hundred bonds per 
[im 2 , which corresponds to area fractions up to 0.01 for 
membrane patch sizes a ~ 5 nm. At these bond area frac- 
tions, the exclusion rejection is negligible. To understand 
the increase in the fluctuation rejection probability with 
decreasing bond fraction X, it is instructive to consider 
how the membrane roughness depends on X. The over- 
all membrane roughness represented by the dashed 
line in Fig. 2(b) diverges with decreasing bond fraction 
X. Wc find that the roughness scales as £_l ~ 0.14&& 
with £f,f, = 1/y/X, see Fig. 3(a). In contrast to the over- 
all membrane roughness, the roughness £x,n at unoccu- 
pied nearest neighbor sites of receptor-ligand bonds (full 
line in Fig. 2(b)) reaches a finite value close to 0.5, or 
Q-ha^/ksT / k in physical units. At small area fractions X, 
the roughness £x,n is proportional to the fluctuation rejec- 
tion probability for hopping steps, see Fig. 3(b). There- 
fore, the membrane rejection probability also reaches a 
finite value in the limit of small bond area fractions. This 
value is about 0.44, which implies that the numerical 
factor by which the diffusion constant of the receptor- 
ligand bonds is reduced due the thermal fluctuations is 
1/(1 — 0.44) ~ 1.8 in the biologically relevant regime of 
small bond concentrations. 

The rescaled bond length z can affect the diffusion if 
it is close to or smaller than the roughness Typ- 
ical receptor-ligand bonds that interconnect cell mem- 
branes have a linear extension l a between 15 and 40 nm. 
For bending rigidities k from 10 to 20 fc^T, this corre- 
sponds to rescaled bonds lengths z — l j 'ay/ n/ksT be- 
tween 10 and 40 at room temperature. For typical bond 
area concentrations X/a 2 in cell adhesion zones between 
10 and several hundred molecules per /urn 2 , the rough- 
ness £j_ — 0.14/y/X is smaller than these rescaled bond 
lengths. The length of the bonds then docs not affect 
the diffusion behavior. However, biomimctic, or 'artifi- 
cial', receptor-ligand systems can be significantly smaller 
than the biological bonds. For example, membrane bonds 
between negatively charged lipids mediated by multiva- 
lent ions such as Cr 3+ can impose local membrane sep- 
arations smaller than 1 nm [15]. Fig. 4 shows how the 
overall membrane roughness and the membrane rejection 
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Fig. 4: (a) Overall membrane roughness £j_ as a function of 
the rescaled length z of the receptor-ligand bonds at three 
different bond area fractions X, and: (b) Fluctuation rejection 
probability Pfi versus z at the same three bond area fractions 
X. For z < 1.3, the roughness £j_ and the fluctuation rejection 
probability Pfi for hopping steps increase with decreasing z . 



probability Pfi for hopping steps depend on the length of 
relatively short bonds. As a function of the bond length, 
the roughness has a minimum at intermediate lengths, 
and increases for small lengths. This increase is a con- 
sequence of the fluctuation-induced repulsion of the mem- 
branes. For small bond lengths z a , the average mem- 
brane separation is then significantly larger than z a , which 
causes an increase in the roughness. This is reminiscent 
of the increase in roughness with decreasing bond length 
in the case of membranes adhering via parallel stripes of 
receptor-ligand bonds [16]. The fluctuation-induced repul- 
sion of the membranes induce also attractive interactions 
between the bonds (see Fig. 5). These entropic attractions 
can lead to lateral phase separation for large bonds that 
occupy several membrane patches [17] or for bonds that 
locally change the rigidity of membrane patches in which 
they are anchored [18], but not for the small bonds consid- 
ered here [17]. Aggregation of receptor-ligand bonds in 
the adhesion zones of cells [19,20] or vesicles [21] has also 
been observed experimentally. In these situations, the ag- 
gregation is presumably driven by a length mismatch (i) 
between bonds and repulsive molecules [21] or (ii) between 




Fig. 5: (Top) Monte Carlo configuration at the concentra- 
tion X — 0.1 of receptor-ligand bonds in the limit of large 
rescaled bond length z . Bonds are shown in black, and the 
local membrane separation is indicated by colors ranging from 
blue (small separations) over green to yellow (large separa- 
tions). The membrane consists of 100 x 100 patches of size 
a 2 . - (Bottom) Monte Carlo configuration at the same concen- 
tration X — 0.1 and the relatively small rescaled bond length 
z = 0.5. The membrane separation z varies between (blue 
color) and 4.85 (red color). The clustering of receptor-ligand 
bonds is caused by the fluctuation-induced repulsion of the 
membranes. 



different types of receptor- ligands bonds [19], which leads 
to a barrier in the effective adhesion potential of the mem- 
branes [18]. 
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Discussion and conclusions. We have focused 
here on the effect of membrane shape fluctuations on 
the diffusion of receptor-ligand bonds between two mem- 
branes. We have found that the fluctuations reduce the 
bond diffusion constant by a factor close to 2 at small bond 
area fractions, compared to planar membranes. Mem- 
brane fluctuations can also affect observed diffusion con- 
stants via the ratio of the actual and the projected mem- 
brane area in the experimental focal plane [22, 23], by 
coupling to the spontaneous curvature of membrane pro- 
teins [22], or by curvature- induced changes of the mem- 
brane thickness [23]. As 'pinning centers' between fluc- 
tuating membranes, the bonds considered here can also 
impede the lateral diffusion of other large membrane pro- 
teins [24]. 

We have assumed that the receptor-ligand bonds in our 
model are rather rigid and have a constant bond length 
l Q . This assumption is justified as long as variations in 
the bond length are significantly smaller than the rough- 
ness £x,„ at bond neighbor patches. At small bond area 
fractions, £± tTl is close to 0.5, which corresponds to vari- 
ations £,_\_, n a,y/kBT / k in the membrane separation up to 
f nm. Another assumption of our model is that Monte 
Carlo diffusion steps of the receptor-ligand bonds imply 
a shift of the local membrane separation and an associ- 
ated change of the membrane bending energy. We have 
previously shown that the Monte Carlo dynamics for the 
membrane separation field z corresponds to a Langevin 
dynamics in the limit of small Monte Carlo step widths 
5z [25]. For the diffusion steps, we have chosen here for 
simplicity a finite Monte Carlo step width that is identical 
with the lattice spacing a of the discrete membrane, which 
prevents the mapping to a continuous Langevin dynamics. 
We have also neglected hydrodynamic interactions, which 
may lead to a further reduction of the diffusion constant. 
Theoretical descriptions of diffusion in membranes mod- 
eled as elastic sheets require assumptions and approxima- 
tions that, in principle, can be elucidated in simulations 
with molecular membrane models [7, 26] . Our theoretical 
results are also accessible to experimental studies of model 
systems with irreversibly bound receptor-ligand complexes 
such as the Biotin-Avidin complex. 
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